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ABSTRACT 

Of the many probes of reionization, the 21 cm Hne and the cosmic microwave background 
(CMB) are among the most effective. We examine how the cross-correlation of the 21 cm 
brightness and the CMB Doppler fluctuations on large angular scales can be used to study this 
epoch. We employ a new model of the growth of large scale fluctuations of the ionized frac- 
tion as reionization proceeds. We take into account the peculiar velocity field of baryons and 
show that its effect on the cross correlation can be interpreted as a mixing of Fourier modes. 
We find that the cross-correlation signal is strongly peaked toward the end of reionization 
and that the sign of the correlation should be positive because of the inhomogeneity inher- 
ent to reionization. The signal peaks at degree scales (£ ^ 100) and comes almost entirely 
from large physical scales (A; ~ 10~^ Mpc). Since many of the foregrounds and noise that 
plague low frequency radio observations will not correlate with CMB measurements, the cross 
correlation might appear to provide a robust diagnostic of the cosmological origin of the 21 
cm radiation around the epoch of reionization. Unfortunately, we show that these signals are 
actually only weakly correlated and that cosmic variance dominates the error budget of any 
attempted detection. We conclude that the detection of a cross-correlation peak at degree-size 
angular scales is unlikely even with ideal experiments. 
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1 INTRODUCTION 

During some epoch between recombination (z ~ 1100) and today, the intergalactic medium (IGM) underwent a transformation from almost 
completely ne utral to almost completely ionized. When and how the IGM was ionized remains one of the most exciting open questions 
in cosmology i Barkana & Loeb 2001; Furlanetto et al. 20061) . The timing and duration of this epoch of reionization contains a wealth of 
information about the first cosmic structures, information which is expected to help explain how the primordial density perturbations observed 
in the cosmic microwave background (CMB) evolved into the complex structures we observe in the low redshift universe today. 

Currently there are only weak observational constraints on the epoch of reionizatio n. The absence of a Gunn -Peterson trough in the 



spectra of 2 < 6 quasars indicates that reionization was complete by z ~ 6 Ipan et al. 2002 



White et i 



al.l2()03l:lFan et al..2006). while integral 



constraints on the total optical depth from WMAP imply reionization began ai z > 10 " jPage et alJl2007l: ISpereel et aklbood). A variety o f 



other methods have been proposed, but none impose strong constraints (for recent summaries, see lFan et alj2006 and Furlanetto et alj2()06l) . 

The 21 cm hyperfine spin-flip transition of HI is the most exciting prospective tracer of the cosmic gas before and during reionization. 
This signal does not r equire the existence of bright background sources, which may be rare at high redshifts, making the entire epoch of 
reionization available /Scott&Reesl ll99d : iMadau et alJll99i iFurlanetto etaPboO ff) . It is a line transition, so that observations at a given 
frequency select out a unique slice of the high z universe. Furthermore, fluctuations in the brightness of the 21 c m signal are caused by d ensity 
(and ionization) inhomogeneities on all scales, making it a direct tracer of the underlying matter distribution 1 Zaldarriaga et aljr2004l) . Thus 
the spectral and angular variations of the 21 cm brightness allow us to reconstruct a 3D map of reionization . The epoch of reionization spans 
the formation of the first luminous sources, which then ionized the surrounding gas jBarkana & LoeblboOlh . Thus a 3D map of the evolution 
of the neutral fraction and density field would provide an unprecedented view of this epoch of structure formation. 

Ionization of the neutral IGM creates free electrons off of which CMB photons may be scattered. This Thomson scattering process has 
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several effects jHu et alJll994l : lDodelson & Jubaslll995h . Photons from multiple lines of sight are blended together, damping the primary 
CMB anisotropies. The scat tered photons also gain some of th e peculiar momentum of the free electrons, generating a secondary anisotropy 
(most recently examined by Giannantonio & CrittendenlboO? ). Finally, the polari zation dependence of Thomson scattering generates a new 
large scale polarization from the anisotropic CMB photon field dZaldarriagall 19970 . 

Taken by itself, the CMB temperature provides only an integral constraint on the column density of ionized electrons, because the net 
damping is nearly independent of the location of the electrons. The induced polarization co ntains somewhat more information about th e 
history of reionization, but it is difficult to extract even in a cosmic variance-limited survey dHolder et alj|2003l : iMortonson & Hull2007n . 
Cross correlating with the 21 cm signal should allow more information to be extracted from both measurements. The cross correlation signal 
arises from the velocity field of ionized baryons, whic h sources the Dopple r anisotropies in the CMB and also traces the linear overdensity of 
neutral hydrogen, the source of the 21 cm brightness. I Alvarez et al. I J2OO6I) showed that combining the two measurements could in principl e 
allow the extraction of the global reionization history, which is otherwise difficult to measure (e.g., Ishaveretalj|l999l : iFurlanettolliooj) . 
One reason is that the signal from the 21 cm radiation will be dominated by foregrounds and detector noise which makes extracting useful 
information difficult. Although many of these same foregrounds appear in the CMB, they are much smaller and more easily removed (and 
of course the thermal noise is uncorrelated). Cross-correlation with the CMB could therefore dramatically improve our confidence in the 
cosmological origin of the 21 cm signal. 

This work is an improved calculation of results presented by Alvarez et al. 1 2006h . In particular, we include the corrections to the 21 cm 
brightness due to redshift space distortions caused by the peculiar velocity field of neutral hydrogen (as in iBharadwaj & Alj2004l) . We employ 
a more physically motivated reionization history together with a new model of the growth of the ionized contrast as reionization proceeds. 
We also re-examine the prospects for detecting the large scale cross-correlation. Unfortunately, we find that it will be extremely difficult 
because of cosmic variance. Throughout this work, we will confine the discussion to large angular scales (~degree) where the details of the 



the ionized bubbles that appear during reionization ( Furlanetto et aU2004a) largely average out. On smaller scales, these arcminute structures 
will also induce interesting cross-correlations (e.g.. 



Cooravll2004 ISalvaterra et al fcooilsiosaretalJboOTh , but the simple models based 



linear theory used here will not suffice to describe them. Also, note that we will assume gaussian fluctuations throughout, as app ropriate 
for lin ear matter fluctuations on the large scales we study. Only on smaller scales will higher-order correlations become important dCoorav 
2OO4I) . 

This paper is organized as follows. In §2 expressions for the 21 cm brightness, CMB Doppler brightness, and the cross-correlation signals 
are derived. Equation ( I22l > is the main result. Section 3 details the reionization model, and §4 presents calculations of the observability of our 
predictions for future radio telescopes. Finally, we conclude in §5. 

Throughout our work the following Fourier convention is used 

3 



/(ft,')) 



/ke 



-jk'n{7jo-77) 



(1) 



The cosmological p arameters are set at f2o = 0.24, Sib = 0.042, JIa = 0.76, h = 0.73, as = 0.8 and n = 0.96, consistent with the 
WMAP year three data 1 Spergel et al. 20061) . and we use the matter power spectrum of Eisenstein & Hij 1 1997 ). We work in natural units 
where c — 1, so that in a flat geometry the conformal time, 



'70 - viz) 



dz' 
Hiz'Y 



equals the comoving distance, r{z) = ryo 
plane wave. 



(2) 

ri{z). Throughout H{z) denotes the Hubble parameter. We will also use the expansion of the 



(3) 



4^^(-i)*j,(fex)y,,„(k)y;„(x) 



(4) 



where Vi is the Legendre polynomial of order I. 



2 21 CM AND CMB DOPPLER TEMPERATURE ANISOTROPY 
2.1 21 cm Temperature Anisotropy 

As long as the spin or excitation temperature, Ts, of the 21 cm transition in a region of the intergalactic medium (IGM) differs from the CMB 
temperature that region will appear in either emission (Ts > Tcmb) or absorption (Ts < Tcmb) when viewed against the CMB. Variations in 
the density of neutral hydrogen would appear as fluctuations in the sky brightness of this transition i Scott & Reeslll99(il : Madau et alj|l997 ; 

Furlanetto et ahlbooeh . 

The optical depth of the IGM in the hyperfine transition is ( lFieIdlll958h 
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T21 = 



3c RAioriff 



-a;Hi, 



(5) 



IQkBvlTsil + z){dvr / dr)' 

where Aiq — 2.85 x lO"'^^ is the spontaneous emission coeffiecient, fca is the Boltzmann constant, uh is the number density of 
hydrogen, and xm is the neutral fraction, dvr/dr is the gradient of the radial velocity along the line of sight, Vr being the physical velocity 
and r the comoving distance. This factor accounts for the bulk motion of the hydrogen, which causes a Doppler shift from the intrinsic line 
frequency, vq = 1420.2 MHz. 

Peculiar velocities induce a small perturbation in the gradient of the proper velocity away from the Hubble flow 



(6) 



where v is the radial peculiar velocity and a{z) is the scale factor. To first order in the density perturbations the optical depth in the hyperfine 
transition is then 



T21 ~ 8.6 X 10"^a;Hi 1 + 5i,{h, j?) 



1 



dVr 



a{z)H(z) dr 



0.02 



/ 0.15 \ / I +z \ 

xn^h^) V 10 J 



1/2 



(7) 



where Tcmb{z) — 2.725(1 + 2) K and we have assumed that the universe is matter-dominated at the redshifts of interest. In arriving at this 
expression we have neglected a number of terms, such as the effect of the CMB dipole in the gas frame and the local effects of gravitational 
redshifting, which also contribute to the 21 cm optical depth. H owever, these additional terms contribute less than 1% to the anisotropy on 
the scales of interest, and equation ^ suffices for our purposes jLewis & ChallinojboOTl ). 

Following Zaldarriaga et al. (2004), we define T2i{n, z) to be the observed brightness temperature increment between this patch and 
the CMB at an observed frequency p corresponding to a redshift 1 -\- z = vq/v ina. direction n, 

^21 (2) ^ T^(r. 



1 + Z 

where we have assumed that T2\ is small. 

We now write the observed differential brightness temperature of the 21 cm emission line at A = 21 cm(l + z) in the direction fi as 



(8) 



r2i(n,2) = ro(z) 



dr?'Vy[r,(z)- 77^21 ( 



n,77 



(9) 



where Vy[r;(z) — rj'\ is a normalized AxW[x\ — 1) spectral response of an instrument which is centered at r;(z) — r;' = 0. T(,{z) is a 
normalization factor given by 



To{z) ~ 23 mK 
and 



0.02 



V'2i(n, 7;) = xm(n,i]) 



0.15 



1 + Sb{n,ri) 



1 + z 
10 



1/2 



dVr 



a{rj)H{vi) dr 



1 - 



Tcmbijf) 



Ts{n,v) 



(10) 



(11) 



We assume that the spectral resolution of the instrument is much smaller than the features of the target signal in redshift space. We thus set 
W[x] = 5°{x), to obtain 



T2i{n, z) = To{z)ij2i[n,ii{z)]. 



(12) 



On very small scales, this is a poor approximation because the finite bandwidth averages over many oscillations, damping the signal. However, 
for our regime of interest at ^ ~ 100 this is expected to be a n excellent approximation (.Zaldarriaga et al...2004.) . 

During most of reionization, we expect Ts 3> Tcmt Ciardi & MadaJ 20031 ; Furlanettol 200€ ). so we neglect the Tcmt/Ts term in 
equation dl lb . By writing the ionized fraction contrast 



(13) 



where Xi is the ionized fraction and the overbar denotes an average quantity, to first order in the density perturbations, equation dl It can be 
written 



i^2i{n,v) = a::Hi 1 + Sb{n,rj) 



dVr 



a{z)H{z) dr 



iiSi{n,ri). 



(14) 



It is important to note here that this is not necessarily a well-defined perturbation expansion. At any given point in space, the ionized fraction 
is either zero or unity, so that averaged over sufficiently small scales (of order the characteristic size of the HII regions). Si itself will be at 
least of order unity. However, we will only examine large scales (^ 10 Mpc) which average over many of these highly ionized regions, so 
equation ( I14t is acceptable. 

Working in Fourier space, the baryon velocity field is related to the density contrast via the continuity equation, Vk = — ik/fc^5i,.k('?), 
where the overdot denotes a derivative with respect to conformal time. At late times on the scales of interest, the perturbations scale as the 
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linear growth fac tor D\, so i5h .k(?7) = Si,,]f^a{'rj) H (rj) f ,wheie f = d In Di/d In a is the dimensionless linear growth rate or redshift space 
distortion factor i Kaisei|[l987 ). Equation l ll4t becomes 



f d'^ 



(27r)3 

The multipoles of the spherical harmonic expansion are then 



aiUz) = "^-^i-iTToiz) 



(27r)3 



{xm{z)Si^{z) Jt[k{rjo - - Xi{z)5a^{z)ji[k{'r)o ~ v{z))]} y/m(k). 



(15) 



(16) 



where jsharadwai & Alj2004l : lBarkana & Loeblboosh 



M^) = [kivo - viz))] + if^^^rrrr^ + ^ IHvo - vim - /7S^4^SfAj'^+2 [^('^o - vim , (i7) 



4£2 + 4^-3 



(2£+ l)(2^ + 3)' 



although we also include the redshift space distortion factor /. This is well approximated by / = f2m^, which for the redshifts of interest is 
approximately one. Neglecting the peculiar velocity pe rturbations amounts t o taking Jiix) — > jtix). 

To this point, our calculation is similar to that of Alvarez et al. i 2006'), except that we have incorporated the velocity correction in a 
different manner. They differentiated the plane wave before expanding in spherical coordinates. This left a factor where fi — n-k/jk| 

is the angle between the line of sight direction n and the photon propagation vector k. On the other hand, we expand the plane wave into a 
spherical basis before taking the derivative. Both approaches are mathematically equivalent, but the advantage of ours will become apparent 
when calculating the cross-correlation. 



2.2 The CMB Doppler Signal 

In Fourier space, the secondary CMB temperature anisotropy from the Doppler effect is given by the line of sight integral 



ro(ri, k) = — Tcmb / drj{—r)fivbi^e 



ikn(ri-rio) ^ 



-Tcmb I d77(-7-)e ^^-^ \e 
lo dv ^ 



d 



ifc(j(i7-r;o)l 



(18) 



We define T{ri) = /™ Av (JTTieQ. SO that r = —aTfiea. Note that with this definition (— r) > 0. 

As above, the baryon velocity field is related to the density contrast via the continuity equation. Thus in real space, assuming an observer 
positioned at the present day (r; = 770), we have 



Expanding in spherical coordinates, the multipole moments are 



afm ^ 4:TvTcmbi-i)'' I dvD(-T)e 
'0 



d^fe 5fc 
(27r)3 W 



^j4fc('7o - v)] 



Yi%i^). 



(19) 



(20) 



2.3 Doppler-21 cm Cross Correlation 

We now calculate the cross correlation power spectrum, C?^~^ 



^af^). We define the 3D power spectrum, {5\i{z)5^^l [z)) = 



(27r)35^(k-k')PM(fc,2)andthe cross-correlation power spectrum between the ionized fraction and density, ((5k(2)'5i.k' iz)} = i2n) S (k~ 
k')Pi»(fc,z). Wehave 



/ 21 / N D*\ 
K'^lmKZjO'trnl 



-Toiz)nmt^^^ I dk[xHiiz)Pssik)Mk{vo - vim 



x^[z)Puik)ji[k{vo - vim I dv'Di{-T)e -^ji[kivo - v') 



(21) 



We perform an integration by parts, and neglect the surface term since at early times r ~ cxa and at late times t{Q), the scattering rate, is 
essentially zero. We obtain 

(aLWaf™) = ~Toiz)nrnb^^^ J^" dri' ^ [Di{-r)e-^] J dkMkii^o - v')] 

X {xmiz)Pssik)Je[k{vo - vim ~ x,{z)P^sik)je[kivo - vim) ■ (22) 

Compared to the expression of lAlvarez et al. we have the factor Jt[k{rjo — vi^))] instead of A/'iji[k{vo — vi^))]- To obtain equa- 

tion l l22t . we have simply required the orthogo nality of the spherical harmonics to evaluate j df2Y«(/i, (f>)Y/,{fi, (f). However, if the factor fi^ 
is included instead, as in Alvarez e t al. 1 20061) . then one must evaluate J dQ,V2ifJ-)Yi{ii, (f>)Yi,{fi, (f>), where V2ip) is the second Legendre 
polynomial. I Alvarez et al.. (2006) approximated this integral with its value for I = Q, introducing the factor 4/3. 
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Ignoring thie velocity corrections to the 21 cm signal for a moment [this amounts to taking Je{x) — > ji{x)}, and using the Limber 
approximation (proved in appendix A): 

k'dk^Mkr)Mkr') - P (fc = ^) (23) 







valid for £^ 1, yields 

2tv ^ ^ 27r dr]^ 



xhi{z)Pss [ -tt] -Xi(z)Pu \ 

r{z) J \r{z) 



(24) 



This form is identical to equation (17) of Alvarez et al. 1 20061) , except without their factor of 4/3. It follows from equation ( I24t that the 



cross correlation power spectrum should roughly trace the shape of the underlying matter power spectrum (at least for uniform reionization, 
or for Si <x 5 SLS below). The matter power spectrum has a broad peak on the scale of the horizon size at matter-radiation equality, keq — 
0.009 Mpc~^ (nmft^/0.128). Using the fact the the conformal time is on the order of 10* Mpc for large z, equation l l24b implies that the 
cross correlation power spectrum will have a peak at degrees scales or ^ ~ fcr ~ 100. For the rest of the paper we take ipeak ~ 100. 

Equation i24i implies another important fact: there is a one-to-one correspondence between k and £, so each multipole samples only 
one scale. In the exact case this is only approximately true; the right hand side of equation ( 123b is not a true delta function but has some 
finite width, which means that modes are mixed. The peculiar velocity corrections further mix the modes as seen through the appearance of 
je+2{kx) and je-2{kx) in equation ( 117b . These redshift space distortions boost the signal by about 20% at ipeak- 

To obtain a numerical result we need a sensible model for the reionization history. On the scales of interest, density perturbations grow 
linearly. We therefore use a growth factor normalized so that Di (zjv) = 1, or 

d / 1 + 2JV ' 



(25) 



Assuming that only hydrogen is ionized, 

nip 

= 0.0525HoQbh{l + zfx,{z), (26) 



where Yp = 0.24 is the helium mass abundance. Taking e ^ « 1, since r ~ 0.1 ( Spergel et al.ll200^ , then gives i Alvarez et al.ll20o3) : 



^ [i)(-f)e-"] = -0.0525m!^i,/i(l + zn)H{z)-^ [x,{z)H{z)] . (27) 



The combination (1 + znY P{k, zn) leaves equations ( 122b and ( 124b independent of the epoch of normalization zn (H^EssH). 



As noted above, the shape of the cross correlation power spectrum for fixed z is determined entirely by the underlying matter power 
spectrum. For fixed I its variation with redshift is determined primarily by the quantity d/drj{DTe^'^), which from equation ( 127b depends 
on d/dz [xi{z)H{z)]. This implies another important fact: the cross-correlation vanishes when d/dz [xi{z)H{z)] is constant and will be 
the largest where [xi{z)H{z)] has the greatest rate of change. Models with faster rates of reionization and reionization occuring at higher 
redshifts will lead to larger cross-correlations (see below). This behavior comes directly from the CMB Doppler signal; unless the ionized 
fraction is changing, the line of sight integral in equation ( |19b suffers severe cancelation due to the oscillatory nature of the Bessel functions. 
Furthermore, this means that if dilution from the Hubble expansion is faster than the reionization rate, then the CMB will be on average 
blueshifted by reionization (and vice versa). During reionization Xi will almost definitely be increasing very rapidly, so redshifting will be 
the dominant effect. 

To get an idea o f how the velocity corrections affect the sig nal we chose a simple parametrization of the neutral fraction to compare 



lo get an idea o t how the velocity corrections attect the sig i 
with previous results i Furlanetto et alj2004bl ; Alvarez et al.ll200a) : 



1 + exp[— (2; — Zr)/Az\ 

Figure [T] shows the cross-correlation at 2 = 15 for homogeneous reionization {Pgi — 0) using equation ( |28b with Az = 1 and Zr = 15. 
The three curves show the Limber approximation of equation ( 123b , the exact result, and the anisotropy ignoring velocity corrections (i.e., 
Ji —> ji), from top to bottom. Note how the Limber approximation works extremely well at £ ^ 100 but overestimates the signal by 
~ 10% near ipeak and fares considerably worse at small £. Velocities b oost the signal on large angular scale by tens of percent but become 
unimportant on small scales. Note that the result of lAlvarez et al. I J2OO6I) corresponds to 4/3 times the lower curve, which at £ 100 is not a 
bad approximation but gets worse as £ increases. 

2.4 The Cross-Correlation Coefficient 

We will find later on that it is convenient to re-cast many of our results in terms of the cross-correlation coefficient r. 



r^/cfCp^, (29) 
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a. 
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Figure 1. Cross-correlation angular power spectrum. The solid curve is calculated using the Limber approximation. The dot-dashed curve is the exact numerical 
integration without the velocity corrections and the dashed curve includes the velocity corrections. We take Zr = 15, Az = 1, and 2 = 15 



where is the angular power spectrum of the 21 cm signal and Cf"'^ is the angular power spectrum of the CMB; in particular, we 
are interested in the Doppler contribution to the latter, Cf^. The maximal cross-correlation coefficient will be rmax ~ \/ /Cf^^ , 
because the 21 cm signal is (nearly) uncorrelated with the primordial CMB anisotropy and other secondaries. 

We can estimate the intrinsic 21 cm signal with a little bit of work fr om equation {1611. Negle cting redshift space distortions, which give 
about a 30% correction at £ — 100, and using the results above, we have I Zaldarriaga et al. 12004 ) 



Cf{z) = {4nfTo{zf J k^'j^je [k{rio - v{z))f [xljP{k, z) + x?P,.(fc, z) - 2xnix,Pu{k, z)], 

where Pa is the power spectrum of 5i . 

We can estimate the CMB power from the Doppler signal from equation J20t . 



D D *\ _ 



TLb dr,|- [b{-i)e-^] j dr;'^ [b{^r)e-^] ^ j dkP{k)Mk{rjo - vmHvo - V% 



(30) 



(31) 



where we have used the orthogonality of the spherical harmonics and integrated by parts twice. Using the Limber approximation (eq. I24t 
we obtain 

'T}e I I p(:;:^] {vo~vf- (32) 



10 / Q 







Note that the amplitude of this expression is primarily determined by the factor d/drj{bfe '^), which as we saw above depends on 
d/dz [xi{z)H{z)]; the Doppler signal will be largest for models of reionization which occur at high redshifts and which proceed quickly. 



3 THE REIONIZATION HISTORY 

The prescription of equation ( 128b is not a well-motivated model of reionization. The quantity of interest for the cross-correlation, the rate at 
which the ionized fraction is changing, is completely determined by the free parameter Az, so the form has no predictive power. Moreover, 
as soon as the ionized fraction reaches 1/2, the rate of reionization begins to slow, so the shape of the signal is prescribed "by hand" 
rather than following from physically-motivated arguments (e.g., recombinations). Our next task is therefore to estimate the signal on more 
physically-satisfying grounds. 

However, in this paper we do not aim to present a complete model of reionization. Instead we will use a simple model that encapsulates 
many of the basic features in the global evolution. We begin by assuming that a galaxy of mass mgai can ionize a mass C'Tigai, where is 
the ionizing efficiency: a measure of the number of ionizing photons produced per baryon inside galaxies. This amounts to the assumption of 
an "inside-out" reionization scenario wherein the gas is first ionized by galaxies within overdense regions. This is consistent with numerical 

. We then associate the reionization rate with the 
star formation rate (iFurlanetta,2006) : 



simulations, at least when the sources resemble star forming galaxies i McOuinn et al J2OO5I 
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Figure 2. Left pane!: The reionizadon model of lFurlanettd j20O6l) with fa 1 at z=l (dot-dashed curve), 10 (dashed curve ) and 15 (solid curve). Right 
panel: The contribution of the homogeneous term to the signal at ipeak for the same models. 



dxi 



c 



d.fc, 



dt 



- aAC{z,Xi)xi{z)ne 



(33) 



where /coii is the mass boun d to halos above some mass t hreshold mmin- In this wo rk we take the minimum mass to correspond to a virial 
temperature of 10^ K tBarkana & Loebl 2001 ) and use the Press & Schechtei 1 1974 ) mass function for simplicity. The second term on the 
right acts as a sink describing recombinations, oa ~ 4.2 x lO"^'^ cm'^s"'^ is the case-A recombination coefficient at 10'' K, C = (ni) /{rie)'^ 
is the clumping factor for ionized gas, and is the average electron density in ionized regions. 

The physics underlying the ionization efficiency is highly uncertain and depends on poorly understood intra-galactic dynamics, so in 
this work we take it as a free parameter. In principle and mmin could vary with redshift due to feedback mechanisms. In most cases, 
feedback decreases and hence the reionization rate. Thus in models that include fe edback, we expect the cross-correlation to be smaller 
and the peak to occur earlier in the ionization history (e.g. Furlanetto & Loebl 2005 ). For example, if the photoheating that accompanies 
reionization suppressed star formation in small galaxies, would decrease significantly once Xi > 0.5, stretching reionization over a longer 
time interval (possibly enough to create an inflection p oint in the history; see below). In e xtreme cases, where recombinations dominate early 
in reionization, ther e could even be a turnover in Xi ([Ce I2OO3I : IWvithe & Loeb l l2003l) . which would strongly affect the cross-correlation 

1 Alvarez etalj2006l) . 

The clumping factor, C, can in principle be computed with numerical simulations that incorporate radiative transfer. However, this 
requires fully resolving small scale structure in the IGM as well as tracking its evolution through reionization. Furthermore, C depends on 
how the I GM was ionized; wheth er low density gas was ionized first or if many photons were consumed ionizing dense blobs when Xi 
was small (Furlan etto & OhI 2005). This means that C must be recomputed for ea ch different set of m odel parameters. Simulations have 
unfortunately not yet reached these goals. We instead use a simple analytic model jMiralda-Escude et al.. 2000) in which voids are ionized 
firstQ This model perhaps underestimates the mean recombination rate, and so reionization probably proceeds a little slower that the model 
predicts. This means that our models will overestimate the signal and will predict a peak that is slightly later. This optimistic assumption 
about the speed of reionization means that our conclusions about the difficulty in observability of the signal should be robust to uncertainties 
in the history of reionization. 

The left panel of Figure[2]shows the ionized fraction as a function of redshift for three scenarios, with reionization finishing by z — 15, 

2 = 10 and z = 7. They have CMB optical depths of r = 0.17, 0.10, and 0.06 respectively. 

We are now in a position to calculate the cross correlation under the assumption that reionization is uniform (Si — everywhere). In 
this case Psi{k) = and only the first term in equation ( 122b contributes. The right panel of Figure|2]shows li^C^^^^ /2Ti for the reionization 
models described above. We show its value as a function of redshift at ^peak; the shape of the angular spectrum is only weakly-dependent 



^ This may at first seem to conflict with our assumption that ionization begins in dense environments. In fact, the ionizi ng photons will begin in larg e-scale 
overdense regions, but within the resulting ionized bubbles they will tend to ion ize the low-den sity gas first. As a result, the lMiralda-Escude et alj jlOOOt) model 
underpredicts the recombination rate by a factor of a few early in reionization jpurlanetto & Oh.2005i) . but that has little effect on our results. 
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on redshift. Each curve corresponds to the matching reionization history on the left. The signal is negative (implying an anti-correlation) 
and peaks where the ionized fraction is changing the fastest. The signal decreases toward lower redshifts for three reasons. First and most 
importantly, the Hubble expansion dilutes the electron densities and hence decreases the scattering rate of CMB photons. Second, /coii//coii 
decreases with cosmic time (though only gently), so the relative ionization rate decreases somewhat in scenarios with late reionization. 
Fina lly, the clumpi ng factor C also increases with cosmic time, slowing the tail end of reionization most severely when Zr is small. 



Alvarez et alj ([200© also used a physically-motivated model for the ionization history that included both galaxy collapse and recombi- 



nations. The major difference with our model is that they prescribed Xi[z) to have a particular functional form that matches smoothly onto 
Xi = 1 (something that we achieve with a dynamically evolving clumping factor). Comparing their Fig. 3 with the z — 15 curve in our 
Figure[2l we predict a larger sign al that peaks later in t he reionization history. This follows directly from the different reionization histories: 
the functional form prescribed by I Alvarez et al. has an inflection point where Xi — 1/2, which causes their estimate to peak slightly 

earlier. As noted above, such a feature could arise in our model from feed back. Our amplitude is also larger because /coii//coii is large 
throughout that epoch; again, the smooth functional form prescribed bv i Alvarez et alj J2006) suppressed this dependence but may be more 
realistic in the presence of strong feedback. 



3.1 Inhomogeneous Reionization 



Reionization is expected to be inhomogeneous for two reasons: the dumpiness of the IGM and the clustering of the discrete ionizing sources. 
The density-ionization cross correlation Psi is a measure of this inhomogeneity. Given that little is known about the epoch of reionization 
and that many of the parameters describing the evolution are completely unconstrained, it is difficult to construct detailed models for this 
correlation. We will therefore rely on the simplest approach that we can. 

To calculate this correlation we need a model of how the perturbations to the ionized contrast Si grow through reionization. Ideally we 
want to relate these to the underlying matter overdensity S so that we can express Psi in terms of Pss- One potential problem is that, on small 
scales, fully ionized bubbles can have Si ^ 1. Fortunately, on large scales these bubbles average out so that linear theory will suffice in the 
regime of interest. 

As before, we assume reionization proceeds "inside- out" so that a galaxy of mass rrigai can ionize a bubble of mass Crrigai around it. It 
follows that, for an isolated large region of overdensity S > Furlanetto et al.lEo04ah . 



^iS)=X,{l + S,)=CfcouiS), 



(34) 



where /coii('5) is the conditional collapsed fraction in our region. Note that the assumption of an isolated region is key here. It is adequate on 
the large, > 100 Mpc s cales we are interested in because ionizing photons are never able to travel such large distances during reionization 
I Furlanetto et alj|2004ar) . A model like that above fails if applied on ~ 10 Mpc scales, because photons from sources outside the region 
become relevant. Expanding this to linear order, we have fco u(S) = /coufl + Here b is the mean galaxy bias (averaged over all galaxies; 
Imo & White^I996) . For an individual galaxy of mass m, the lPress & Schechten ( 119740 mass function yields 



6(m) 



1 + 



cr(m)2Di(z) Sc{z)Di{z) 



(35) 



where 5c is the critical density for spherical collapse and (m) is the fractional variance of the density field when smoothed on spheres of 
mass m. The average over all galaxies is 



b = — !— / dm ^n{m)b{m) 

./coll J P 



(36) 



Furlanetto et al 



2004ah by integrating over the galaxy population 



Note that the same expression follows from the reionization model of |] 
of each ionized bubble in a large-scale region of overdensity S. We show the mean bias as a function of redshift in the upper left panel of 
Figure [51 

In this model, we therefore have Si — bS. Before proceeding further, we must address two more potential problems with our approach. 
The condition that the ionized fraction everywhere be less than unity implies 



1 > X, 



Xi{l + bS), 



(37) 



which can be written 



(38) 

Xi 

On the scales of interest we can estimate 5 ~ a{R). I — 100 corresponds to ii ~ 200 Mpc, where (t(200 Mpc, z — 10) ~ 0.004. We are 
interested in the epoch in which Xi ~ 0.5; so b can be as large as ~ 360 and our approximations will still remain valid. 

This prescription poses one final problem: as reionization proceeds, the ionized contrast grows monotonically with the bias, irrespective 
of the ionized fraction. This ignores two subtleties of reionization. First, even on these large scales the ionized features will eventually grow 
large enough that the regions are no longer independent (i.e., low-density voids will eventually be ionized by their neighbors). Moreover, 
iFurlanetto & Oh C2005i) have shown that recombinations become important for larger bubbles. They imprint a maximum size -Rmax on 
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reionized regions; any extra photons produced in thiese bubbles are canceled out by recombinations. Although this process is difficult to 
model in detail (especially given the weak constraints on structure in the IGM at high redshift), we can approximately account for these 
effects by modifying Si to 



5,; = 6<5(l-/sat), 



(39) 



where /sat is the fraction of the volume in the "saturated" bubbles (with R > -Rmax)- As Xi approaches unity, so must /sat, so the effective 
bias 6(1 — /sat) eventually goes to zero. In other words, we subtract off photons produced inside regions where recombinations dominate, 
because they are consumed before helping to ionize new material. Unfortunately -Rmax depends sensitively on the distribution of high-redshift 
Lym an-limit systems, whose origins are uncertain. There is thus no particularly good model for -Rmax or /sat • However, the curves in Fig. 
10 of lFurlanetto & Olj jlOOSh are reasonably well-described by 



/sat(a;i) = exp 



(40) 



0.2, 0.1, and 0.05. These values of 
= 9, and z — 12 respectively (with 



where Axi is a parameter which encodes the effects of the recombining bubbles (which increases as i?max decreases, because the saturation 
limit is reached earlier for smaller bubbles). T he upper right panel of F ig. [3] shows /sat for Axi 
Axi are chosen to roughly match the results of Furlanetto & OhI i 2005 ) and correspond to 2: = 6, 
recombinations setting in earlier at lower redshifts because the dumpiness C increases with time). 
Finally, the density-ionization cross correlation power spectrum is given by: 

= b{z)il~ Ut)Pik,z). (41) 

Current simulations do not have large enough box sizes to test this approximation on the > 100 Mpc scales of interest. However, on the 
largest scales that we can compare (A: ~ 0.1 Mpc~^) , our result is ~ 10% larger than that obtained from "semi-numerical" simulations of 
the reionization process, at least early in reionization jMesinger & Furlanettdl2007 ). As reionization proceeds, the simulations show that the 



amplitude on this scale gets even larger; this is because the characteristic size of ionized bubbles reaches k ^ 0.1 Mpc^^ about midway 
through the process. We do not expect such an amplification to occur at the scales relevant to the cross-correlation so long as the bubbles do 
not reach such large sizes. 

The lower panels of Figure [3] show the variation in redshift space of the cross-correlation signal, with the left panel showing the 
individual contributions of the homogeneous and inhomogeneous terms and the right panel showing the sum. For reference, 2; = 7, 10 and 
15 correspond to 1/ =177.5, 129.1 and 88.8 MHz respectively. 

The bottom left panel of Figure [3] shows that the positively correlated inhomogeneous part of the signal strongly dominates over the 
anti-correlated homogeneous part. This is because the inhomogeneous component is so strongly biased (upper left panel). Its amplitude 
increases rapidly with redshift for two reasons: the quantity d/dz[H{z)xi{z)] is larger for earlier reionizati on and the galaxy bi as grows 
large at higher redshifts (see the upper left panel of Fig.[3}. Our amplitudes are somewhat larger than those of lAlvarez et al.l ( l2006h because 
our reionization models allow larger dxi/dz and because our cosmology yields larger bias factors at high redshift. 

In the lower right panel we explore the effects of varying the rate at which bubbles saturate by varying the parameter Axi. If saturation 
sets in earlier the peak amplitude of the signal decreases, and the peak also moves earlier in reionization. However, in all of these cases /sat 
has little effect except in a brief redshift interval. This is simply because Xi increases from ~ 0.5 to ~ 1 so quickly. We find, therefore, 
that measuring the peak amplitude of the cross correlation does not precisely yield the rate at which reionization is occurring (or, in our 
simple model, roughly the rate at which gas collapses onto galaxies); rather, it provides the rate at which photons are able to escape from the 
environs of their source galaxy. Late in reionization, the overall amplitude will be suppressed by the recombining bubbles that host most of 
the ionizing sources. Recovering Xi{z) will require carefully modeling this effect. 

Wh en inhomogeneous r eionization is included, our models robustly predict a positive correlation between the CMB and 21 cm signals 
(see also Alvarez et al.l2006l) . This follows from the strong clustering of the first galaxies and reflects our assumption of "inside-out" reion- 
ization. An anti-correlated signal would imply that low-density gas is ionized at least as quickly as high-density gas; thus the overall sign of 
the correlation would provide a clean test of these contrasting approaches to reionization. 



4 FOREGROUNDS AND OBSERVABILITY 

There are two unavoidable sources of noise in any attempt to measure cosmic signals: the intrinsic noise associated with the detector and the 
error associated with sample variance (especially at low multipoles). The latter results from having only one sky and only 2£ + l independent 
samples at each £: 

where /sky is the fraction of sky being observed. For a single 21 cm observing field this is given by fsky ~ dVl/A-K ~ }? / {AdA-K), 
approximated by the diffraction limit of the telescope, where Ad is the area of a single dish and A = 21(1 + z) cm. 
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Figure 3. Upper Left: The mean galaxy bias function. Upper Right: /sat for Ax; = 0.05 (solid curve), 0.1 (dashed curve), and 0.2 (dot-dashed curve). 
Bottom left: The homogeneous and inhomogeneous contributions to the cross correlation at Ip^ak corresponding to reionization models with Xi=i at z as 7, 
10 and 15. In each case, the positive (negative) curves show the contribution from Pgi (P). All curves assume an /sat with Axi = 0.05. The corresponding 
reionization histories are shown in Fig.|2] Bottom right: The total signal at ip^ak- At each z, we show results for the different values of the parameter Axi in 
/sat shown above. 



While the cosmic variance contribution is purely geometric, the detector noise depends on the details of the observation. For 21 cm 
experiments, the RM S detector noise fluctuation per visibility of an antennae pair observing for a time to in one frequency channel is 
i McOuinnetallbood) 



AV"" = 1!^ , (43) 

where Tsys{i^) is the total system temperature. An irreducible limit on Tsys is provided by the radio sky itself, so Tsys > Tsky In fact, at 
the low frequencies relevant to the redshifted 21 cm sky, the noise is strongly dominated by synchrotron electrons from fast electrons in the 
Milky Way, so to a good approximation Tsys ~ Tsky A good rule of thumb for the sky temperature of high latitude, "quiet" portions of the 
sky is (iFurlanetto et al..2006i) 
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Figure 4. The thermal detector noise for 21 cm observations with the MWA (dot-dashed curve) and the SKA (dashed curve) and the intrinsic 21 cm signal for 
the reionization models in Fig.|2](solid curves). 



'sky 



180 



V 180MHz J 



It follows that the detector noise in the auto-correlation for a single baseline is 



Ad J Auto 
The time spent observing each visibility is 

Adto 



A2 



-n(u) 



(44) 



(45) 



(46) 



where u = £/2tt and n(u) is the average number of baselines that can observe the mode u a t any instant, normalized so that its integral over 
the u plane equals the total number of baselines in the interferometer i McOuirm et alJllOO^ . The covariance matrix for an interferometer is 
then 

2 



c 



21. N 



(u) 



Ad J AvU 

The total error on the autocorrelation is a combination of sample variance and detector noise 

3 2 1 ^ ^^^"^ 

A_\ ^ays , /'^21\2 

Ad) Autoniu) ^ ' ' 



ACT 



(2^+1) /sky 



(47) 



(48) 



We estimate using equation ( 1301 ) and the inhomogeneous reionization model described above; note that Pii{k) = P(l — /sat)^f'(fe) 
on these large scales. Figure |4] shows the intrinsic 21 cm signal for our fiducial reionization models as well as the thermal detector noise, 
£2^21, jv j-^j. j.^^ planned experiments, the Mileura Widefield Array-Low Frequency Demonstrator (MWA) and the Square Kilometer Array 
(SKA). We have assumed a frequency window of Au = 1 MHz, 1000 hours of observing time, and I = 100. The area of a single dish is 
taken to be the minimum baseline: (4 m)^ for MWA and (16 m)^ for the SKA, and we assume a filled core outside of which the density of 
antennae falls like r~'^ . The total effective areas are 1 km^ for the SKA and ~ 7 x 10~^ km^ for the MWA. 

The 21 cm signal falls sharply early in reionization because of the prefactor on the power spectrum, [2;hi(z) — b{l — fsB.t)xi{z)f' . 
These two terms cancel at a point early in reionization when the increased total density matches the decreased ionized fraction; at later times, 
the inhomogeneous term dominates and the signal increases again. Obviously, except when reionization occurs at 2 = 15 the intrinsic 21 
cm signal significantly dominates the thermal noise. Thus any measurement of the cross-correlation is sample variance limited. Note the 
difference with typical 21 cm experiments, which are usually dominated by thermal noise: in this case, we are interested in such large scales 
that the 21 cm telescopes can measure the modes in detail at least neglecting foreground contamination). We neglect the 21 cm detector noise 
for the remainder of the paper. 

For the CMB, experiments like WMAP are already nearly signal-variance limited on scales near Ip^ak- Thus we neglect the thermal 
noise term entirely for the CMB. 
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One of the general advantages of cross correlation is that any uncorrelated noise (such as the thermal noise in each detector) will vanish. 
Unfortunately, many of the strong foregrounds for 21 cm observations also appear at CMB frequencies (such as the Galactic synchrotron 
radiation and r adio point sources). However, we will assuine op t imistically that cle a ning of the separate maps will effectively remove this 
contamination dZaldarriaga et al.ll2004l : [Morales & HewittlEoM : ISantos et alj|2003l : [Morales et alj|2006l : iMcOuinn et alJbood) : in reality, 
foregrounds could make the signal even more difficult to detect than estimated here. 

We can then estimate the 1-a error for the cross-power spectrum from 

1 



{21 + l)/.kyA^ 



where Al is the size of the bins over which the power spectrum is averaged (l - 
coefficient, the expected fractional error (eq.|49t can be writtenlf] 



c 



(49) 



< I < (, + M/2). In terms of the cross-correlation 



(50) 



\/(2^ + l)/skyA^ r 

This result has a simple physical interpretation. The first factor on the right hand side is the fractional error from cosmic variance for any 
autocorrelation. If the two signals are weakly correlated, r ^ 1, only a portion of each intrinsic signal is useful for the measurement and the 
cosmic variance limit increases accordingly. This can be recast as a constraint on the minimum correlation coefficient r observable with a 
given signal-to-noise on the cross-power: 

-1/2 



r > 



/ACf 



21-D 
I 



(2£+l)/ekyA^-l 



(51) 



Thus the observability of our signal will depend on the cross-correlation coefficient. We can estimate its value by comparing equa- 
tion ( 132b for our reionization models to the total primordial CMB anisotropies. For the models in Figure|2l we find l?Ci"^ /2tv = 19.8 /.tK^, 



2.4 ijK\ and 0.4 ^jK" for the z = 15, 



10, and z = 7 cases respectively. The primary anisotropy is « (40 /xK) at £ = 100, so even 



if the Doppler component of the CMB is completely correlated with the 21 cm signal (see below), we must still have r < 0.11, 0.04 and 
0.015 for the z — 15, 10 and 7 cases, respectively. One can already see that the prospects for detecting the cross-correlation are dim: the 
CMB at ^ ~ 100 is well-fi t by the primary ani sotropies alone, ignoring the Dopper terms from reionization, and has been observed to nearly 
the cosmic variance limit I Spergel et al.ll20o3> . The fractional cosmic variance errors for C^^'^ will be larger than those on the CMB by a 
factor 1/r. 

For a concrete (though optimistic) estimate, let us assume that we have /sky = 1/2 and sum over all multipoles. Thus we forsake 
any shape information in order to test the null hypothesis, that is, whether or not correlation can be distinguished from no correlation. 
For this calculation, we assume for simplicity that the correlation coefficient between the Doppler anisotropy and the 21 cm brightness, 
ry^C^^^^ /Cf^~^ , is independent of scale and we choose the obse rved frequency to correspond to the peak of the cross- 



r2i-D 



correlation. We can then find the total signal to noise ratio from (compare eq. 5 1 of Slosar et alj2007 ) 



'21-D 
I 



2 

r2i-D 



^(2£+l)/,ky 



^cmb ' 



(52) 



since Cf^ /C1"^^ is small. We find signal to noise ratios of SNR ^ 2 r2i~D, 5 r2i-o, and 14 r2i-D for reionization ending at z — 7, 10, 
and 15, respectively. If the 21 cm signal and the Doppler signal were perfectly correlated, then we could conceivably detect the correlation in 
all these scenarios. However, in reality the correlation coefficient between the 21 cm signal and the Doppler signal is r2i-D ~ 0.3. Thus it 
is unlikely that we can detect the cross-correlation if reionization ends at 2 < 12. For reionization at higher redshifts, a detection is possible 
provided we can view > 50% of the sky. 

Of course, the above estimate assumes a single frequency band for the 21 cm map; increasing the frequency coverage will provide 
more samples and improve the total signal-to-noise (though at the cost of losing any information about the signal's evolution with redshift). 
However, i ~ 100 corresponds to 200 Mpc. Demanding that each 21 cm frequency slice be statistically independent then requires that 
the slices be ~ 10 MHz thick, which corresponds to Az ~ 1. Thus there are in fact relatively few independent slices available, unless 
reionization is considerably more extended than in our simple models (in which case each individual slice would have a smaller signal 
anyway). 



5 DISCUSSION AND CONCLUSION 



In this work we have calcul ated the large scale c ross correlation between the CMB temperature anisotropy and the 21 cm background, 
improving the calculation of Alvarez et al. 20061) . The signal is expected to come from large angular scales (£ ~ 100) corresponding to 



Note that our definition of the cross-correlation coefficient r ignores the noise terms. 
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~ 200 Mpc at redshifts z ~ 10. On these scales hnear theory is still valid, which greatly simplifies the analysis. We have also presented 
a new model for the evolution of the ionized contrast as reionization proceeds that includes the effects of recombinations in limiting the 
apparent bias of reionized regions. 

The cross correlation arises from the connection between linear overdensities and the baryon velocity field via the continuity equation in 
linear theory. Fluctuations in the underlying density field source brightness fluctuations in the 21 cm background, while the baryon velocity 
field Doppler scatters the CMB. The evolving ionized fraction, which is analogous to the evolving gravitational potentials in the integrated 
Sachs-Wolfe effect, reduces the line of sight cancelation that is usually associated with Doppler contributions to secondary anisotropics. 
Moreover, the redshift information contained in the 21 cm signal (as a binning in frequency) allows us to reconstruct 3D information about 
the cross correlation. 

Our calculations include improvements to the approach originally presented by Alvarez et al. ( 2006h . The key differences are the treat- 



men t of the bu lk velocity corrections to the 21 cm signal, the reionization model, and the ionized fraction-density cross-correlation. 



Alvarez et al. (2fl06) chose a particular functional form ln[l — Xi] — —C^o{z)fcou{z) for the reionization history. Although it is driven 
by the same mechanism as ours early in reionization (the time evolution of /coii), this form demands that the reionization history slow when 
Xi > 0.5. Although such a delay can be loosely attributed to either feedback or recombinations, there is in fact no requirement that they set 
in so early, or that they slow the subsequent evolution so dramatically. 

Our model, on the other hand, explicitly relates the collapsed fraction t o the ionized frac t ion, to which we add the effects of recombi- 
nations (which limit the bubble sizes and mediate the source bias). Following Furlanetto & OhI 1 200 5*) . we assume that ionized bubbles grow 
until they reach a size i?max, where recombinations prevent subsequent growth. Any more ionizing photons produced in these bubbles are 
canceled out by recombinations. As reionization proceeds more and more bubbles saturate and no longer contribute to the global evolution 
of the ionized fraction. As this process unfolds, more of the collapsed matter enters these saturated bubbles, so the ionized contrast between 
over- and underdense regions begins to get smaller, eventually shrinking to zero as all regions become ioniz ed. 

These differences change the amplitude and evolution of the signal but not its qualitative behavior. Like Alvarez et al. i 2006h . in the case 
of homogeneous reionization we find an anti-correlation between the 21 cm signal and the CMB. In this case an overdensity simultaneously 
increases the brightness of the 21 cm background and cools the CMB through Thomson scattering across regions where Xi is increasing. The 
contribution from the homogeneous term decreases with redshift due to dilution of the matter density through Hubble expansion. 

The result is markedly different when inhomogeneous reionization is included: the 21 cm and CMB signals become strongly correlated. 
Ionized clumps still scatter and cool the CMB photons through the Xi gradient. However, these regions also cause a significant negative 
fluctuation in the 21 cm brightness due to the large deficit of neutral hydrogen. This term is strongly boosted by the galaxy bias and dominates 
over the anti-correlated homogeneous reionization contribution. 

Of course, this positive correlation rests on our assumption that reionization proceeds inside-out. If instead low density gas is ion- 
ized first, the inhomogeneous component would become anti-correlated: the underdense regions would blueshift the CMB photons while 
decreasing the 21 cm brightness both through a decreased matter density and ionized fraction. 

Unfortunat ely, we have shown a bove that the 21 cm and CMB signals are too weakly correlated for the cross-correlation to be easily 
detectable. Like Alvarez et alj ( 20061) . we find that thermal noise will not pose a significant problem because the correlation only appears 



on large scales, where 21 cm telescopes are sample variance-limited. However, when we include cosmic variance of the 21 cm signal, the 
prospects for detection look dim. The cross-correlation coefficient between the total CMB signal and the 21 cm background is only r < 0.03: 
it is small because the Doppler contribution makes up < 1% of the total CMB anisotropy, and only ^ 30% of that correlates with the 21 cm 
background. Experiments targeting this signal will require extremely large fields of view, and even then will be unable to recover the signal at 
high signal-to-noise or to trace its angular power spectrum. Prospects are actually best for early reionization, because then the CMB Doppler 
term is increased by the denser gas during reionization. (Note also that we have ignored foreground contamination in the 21 cm signal, which 
must be cleaned and will further degrade the detectability.) 

Although it is possible that our models under-predict the cross-correlation signal, our predictions for the detectability is robust. We 
have shown that the signal to noise ratio is limited primarily by the smallness of the Doppler signal. The size of the Doppler contributions 
is determined entirely by the rate of reionization and the epoch at which reionization occurs. Our models already provide relatively fast 
reionization, so the signal can probably only be increased in this fashion by a factor of ~ 2, which does not change our conclusions. 
Furthermore, while detection of the cross-correlation at £ ~ 100 is limited by the cosmic variance of the 21 cm and CMB signals, it does 
not help to move to higher multipoles where there are more modes available. On the one hand, moving to higher I may increase the effective 
bias, b{z), and thus the cross-correlation, as the scale approaches the characteristic bubble size. However, this aspect would only boost the 
21 cm signal (which does not hel p much because it also increases the cosmic variance) and not the Doppler signal. In fact, as shown by 
Giarmantonio & CrittendenI 2007 ) the Doppler signal, and hence the ratio C^^ /Cf^^^ , begins to sharply decrease after I ~ 100: any 



gains made in sampling more modes are immediately lost. 

We thank M. Zaldarriaga and M. Alvarez for helpful comments on the manuscript and A. Mesinger for sharing simulation data with us. 
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APPENDIX A: INTEGRAL APPROXIMATION 

Consider the integral representation of the 3D Dirac delta function, 

(Al) 

Expanding the plane wave in spherical coordinates, we obtain 

"' ' " ^ (A2) 



(27r)^(5(r - r') = j d^fcexp[ik ■ (r - r') 



(2^)^5(r-r') = {^Tvf j d='fc^(-i)*j,(fcr)y,„(k)y/„(i^) ^(i)^'j,,(fcrOy^w(k)yr™'(r') 

£' ,m' 

= (4^)' 51 f k'dkje{kr)ji{kr')Y,^{r)Y;^{r'), (A3) 

where we have used the orthogonahty of the spherical harmonics. Now in spherical coordinates we can express the Dirac delta function as 
S{r~r')=S{Q,n') ^^'^^/\ (A4) 
so that multiplying equation iA3l by Y{m(f ') and integrating over dQ' yields 
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^ k'dkMkr)Mkr') = ^^-^^-^. (A5) 

In this paper we encounter an integral similar to the left hand side of equation l |A5l l, but with another smooth function of k, P{k) /k^,m the 
integrand. The spherical Bessel functions je{kx) are very small for kx < £ and start to oscillate for kx ^ I. As can be seen from l lASI l, the 
Bessel functions are out of phase for any separation of the points r and r' . The integral will thus only receiv e contributions from a r egion 
around the first peak, which occurs at ~ i/x. We then make the approximation P{k) ~ P(k — £/r) (e.g.. lZaldarTiaga et al ] |2004h . and 
pull it out of the integral to obtain the approximation used in the main body of the paper: 



,2 „ P(k) „ /, £\ 5(r-r') 

k^<ik^ji{kr)ji{kr') P (^fc = -j ^ ' . 



(A6) 
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